Edge states interferometry and spin rotations in zigzag graphene nanoribbons 
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An interesting property of zigzag graphene nanoribbons is the presence of edge states which are extended 
along its borders but localized in the transverse direction. We show that because of this property, electron 
transport through an externally induced potential well displays two-paths-interference oscillations when subject 
either to a magnetic or a transverse electric field. This effect does not require the existence of an actual 'hole' 
in the nanoribbon's geometry. Moreover, since edge states are spin polarized, having opposite polarization on 
opposite sides, such interference effect can be used to rotate the spin of the incident carriers in a controlled way. 

PACS numbers: 72.25.-b,73.23.-b,85.75.-d,81.05.Uw 



Graphene, a two dimensional array of carbon atoms in a 
honeycomb lattice, is a veryinteresting material with unusual 
electronic properties 111 |2l [sl |4|] . It has attracted much of at- 
tention since its first experimental realization as it offers 
a great potential for technological applications while, at the 
same time, it has lead to the observation of new physical phe- 
nomena such as an anomalous quantization of the Hall effect 
i?, 8], observable at room temperature, or the manifestation of 
the Klein tunneling paradox in transport iHlillTol], among oth- 
ers. The key for understanding graphene's peculiarities relies 
on its band structure: electronic excitations around the Fermi 
energy (Ep) can be described by an effective Hamiltonian that 
mimics the Dirac equation for massless chiral fermions where 
the spin is replaced by a pseudospin (the two inequivalent sites 
of the hone yco mb lattice) and the speed of light by the Fermi 
velocity ||4ill ll] . The actual spin plays no crucial role in bulk 
samples. 

A novel effect unique to graphene appears in graphene 
nanoribbons (GNRs): when the termination of the GNR cor- 
responds to a zigzag ordering of the carbon atoms (see Fig. 
la) the system presents edge states [12., .isi [l4ll . That is, 
there are eigenfunctions that are extended along the zigzag 
nanoribbon (ZGNR), but that decay exponentially away from 
the edges towards the center of the ZGNR. These estates have 
recently been observed in graphite surfaces near monoatomic 
step edges [15]. From the theoretical point of view, they can 
be easily obtained from either a discrete tight-binding model 
for the honeycomb lattice lITH ITill or a low energy effective 
Hamiltonian (Dirac equation) fl'V]. If only nearest neighbors 
hopping is considered in the former, the edge states have an 
exponentially small group velocity Vg, which leads to a high 
density of states near the Ep of the undoped material. These 
states have been studied in detail by several authors (see [4] 
and refs. therein) including the recent proposal of a novel 
quantum spin Hall effect in the presence of spin-orbit coupling 
1116(1 . When next-to-nearest neighbors hopping is taken into 
account, the edge states become dispersive — they acquire a fi- 
nite Vg — and more stable fl7]. In addition, electron-electron 
interactions lead to a magnetic ordering of the edge states 
lll2i iTsIl and the appearance of an energy gap in the band 
structure ifisll . Since the resulting edge states are then spin- 



polarized, several groups have proposed to use them for spin- 
tronics applications such as creating pure spin currents [19,1 or 
inducing half-metallic behavior with electric fields lIlSll . 

Here, we analyze electron transport through a ZGNR with 
a potential well (PW) created by external gates and tuned in 
such a way that transport inside the well is governed only by 
the edge states. In this case, while the current flow is essen- 
tially homogeneous outside the PW region, it flows along the 
edges inside it. We show then that the system behaves as 
a two-paths interferometer, even though the ZGNR is struc- 
turally homogeneous, an effect unique to the ZGNR band 
structure. Interference between the two paths can be tested 
by either using a magnetic or a transverse electric field to tune 
the orbital phase difference between the two branches. 

Furthermore, since the ground state corresponds to an an- 
tiferromagnetic ordering of the polarization of the two edges, 
each path corresponds to a different spin orientation. Then, if 
the initial spin polarization of the incoming electron, set for 
instance by a ferromagnetic contact, is perpendicular to the 
intrinsic polarization of the ribbon, the two-paths interference 
leads to a rotation of the spin of the carriers whose angle can 
be controlled externally, offering an interesting potential for 
spintronics. 

We describe the ZGNR in the tight-binding approximation. 
The Hamiltonian then reads H = ^^qnr +^cxt + ^^int, where 
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describes the ribbon. Here, (6|^) creates an electron on a 
Wannier orbital centered at site ri of the sublattice A (B) with 
spin (T, i ~ 2.8eV and t' ~ — O.lt f4f] are the nearest- and 
next-to-nearest- neighbors hopping parameters, respectively. 
The symbols (...) and ((...)) restrict the sum to the corre- 
sponding neighboring sites. The borders contain A sites on 
one edge and B sites on the other. iJcxt, which describes 
the action of external gates, is defined below. Finally, tl\^t 
describes the electron-electron interaction. Because of the 
high density of states induced by the edge states, the system 
is magnetically unstable. DFT and Hartree-Fock calculations 
I la |20(] show that the ground state coiTesponds to an antifer- 
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FIG. 1: (a) Scheme of a ZGNR. The energy of a 32-ZGNR as a 
function of the wavevector along the x axis is shown for: (b) Ba ~ 0; 
(c) Ba = t'/2; (d) Ba = t', (e) Ba = 1.2t'. The bands connecting 
the two non-equivalent Dirac points correspond to the edge states. 



romagnetic ordering of the sublattices' magnetization. Since 
the latter is mainly localized at the edges, and to capture the 
essence of this effect, we take into account such interaction by 
introducing an effective magnetic field only at the edges sites. 
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where a {(5) labels the top (bottom) edge. We take this field 
to be perpendicular to the plane of the ZGNR (z axis). In 
the ground state the two edges have opposite magnetizations, 
— —Ba- The value of Ba should, in principle, be deter- 
mined by a self-consistent calculation. Since its precise value 
depends on the chemical passivation of the edges 112 U 12211 . and 
in order to discuss different situations, we take it here as a free 
parameter |23]. 

Figure [T] shows the energy dispersion of a 32-ZGNR ll2411 
for different values of Ba- Several bands originated from 
the quantization along the y axis are clearly visible. The 
bands in the range k^a G [27r/3, 47r/3] that are close to the 
Dirac point, E « 3i', are the ones that coiTespond to the 
edge states with a characteristic localization length \{kx) ~ 
-3ao/21n|2cos(fc^a/2)||[l3]. Here, a = y/Sao is the lat- 
tice parameter with ao the C-C bond length. For Ba = 
(Fig. [TJ)), there are two of those bands (for each spin ori- 
entation) that are almost degenerated — there is an exponen- 
tially small gap between them. They essentially correspond 
to the symmetric and antisymmetric combination of the ex- 
ponentially decaying solutions of each individual edge. For 
Ba ^ (Figs. [Itmillt), both the spatial and the spin de- 
generacies are broken. For each spin orientation, each band 
now corresponds to states localized on a different edge. The 
energy dispersion is approximately given by E{kx) — it' + 
[t' ± jjL bB a) {"2, cos kxa + 1) — note that it is nonzero due to 
the nonzero value of either t' or Ba |I7, 25tl. The key point is 
to notice that, for a given energy, the states with opposite spin 
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FIG. 2: (Color online) (a) Schematics of the proposed setup and the 
potential profile (which depends only on x). The two-paths character 
of the current flow allows for interference effects to manifest; (b) 
Conductance of a 32-ZGNR (W = 45ao) as a function of Vg for 
different values of B^ and Ev = -0.6f, L = 2400a, A = 30a 
and Ba ~ Bb — 0- G only changes in the presence of edge states. 
Inset: G as a function of B± for Vg/t = -0.384 (■), -0.356 (A) 
and -0.31 (♦). 



polarization in the z direction are localized on opposite edges 
oftheZGNR- 

Let us now consider the transport properties of a ZGNR in 
the presence of an electrostatic potential created by external 
gates. 



(3) 



where f{x) is a smooth function describing a PW of height Vg 
(see Fig. |2^). For simplicity, we use a sum of Fermi functions, 
with the parameter A playing the role of the temperature, to 
set the spatial profile of f{x) (which depends only on x)- We 
assume that Ep < 3t' < far from the PW which ensures that 
the cuiTent caiTying states in that region are extended through- 
out the entire width of the ribbon. On the other hand, Vg < 
can be tuned in such a way that Ep ~ Vg corresponds to the 
energy of an edge state. For the sake of simplicity, we dis- 
cuss first the conceptually simpler case Ba — 12M- Then, if 
f{x) changes smoothly , the electrons' wavefunction will adi- 
abatically change from extended to localized, while keeping 
its band index and having a position dependent wavevector 
kx(x)- CoiTespondently, the charge flow will 'split' in two 
paths inside the well and merge again afterwards, creating a 
'hole' in its spatial distribution (Fig. |2^). In this way, we have 
create an interferometer which can be tested by introducing a 
relative phase difference between the two paths. 

As the Aharonov-Bohm effect provides the simplest way to 
do this, we introduce a magnetic field B± perpendicular to 
the ZGNR (via a Peirls substitution in the hoppings) and cal- 
culate the zero temperature conductance using the Landauer 
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approach fz^. For that, we separate the system into a cen- 
tral region (containing the PW) and the lead regions and use 
the standard recursive method to obtain the lattice Green func- 
tions Ii26il27fl and the transmission coefficient from them. Fig- 
ure |2j) shows the conductance G of a 32-ZGNR as a func- 
tion of the Vg for different values of B±. It is apparent that 
G changes with B± only when Vg is below the threshold 
where the edge states participate on transport (indicated by 
the arrow). The inset shows the oscillatory behavior of G 
as function of B± for three different values of Vg. The pe- 
riod is roughly (po/A' ~ 1.3T with the flux quantum and 
A' ~ LcsW with Lcs ~ (L — 4 x 3.5A). An increment of 
the period, due to the reduction of the effective 'hole' area, is 
difficult to see since the visibility of the oscillations is rapidly 
lost. In addition, and despite this seemly simple picture, the 
behavior of the conductance is more involved as it shows pro- 
nounced narrow dips when B± ^ 0. This is related to the fact 
that bonding and antibonding bands are mixed by (recall 
that for i?a = the gap is exponentially small) and then both 
bands get involved in transport which in turns leads to Fano- 
like interference between them 02811 . 

A more interesting situation occurs for Ba ^ 0. As we 
mentioned above, in this case, both the spatial and the spin 
degeneracy are broken. Therefore, an incoming electron with 
its spin quantize along the z axis, will follow either the up- 
per or lower path (colored arrows in fig. |2^) depending on 
whether its spin is 'up' or 'down'. Clearly, in this case there 
is no interference and the conductance is independent of B^. 
Nevertheless, it can be readily verified that if the incoming 
electron is polarized in the fi = cos tfx + sin y direction — 
its spin state being denoted by \n "f ) — it will be rotated 
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where ^ is the relative phase of the transmission amplitude of 
the two paths. Due to the symmetry of the setup, the spin pro- 
jection remains on the plane of the ZGNR. The probability for 
an electron to keep its spin orientation is cos^(^/2) and so we 
expect the conductance between two collinear ferromagnetic 
leads |29] to oscillate as a function of ^. Note that we have 
assumed that L <C icorr, where icorr is the spin correlation 
length of the ferromagnetic order along each edge llsoll . 
Figure [3] shows the spin-resolved transmission probability 



Ta+ for an incident electron with spin 
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transmitted with spin cr = ± (in the same axis) as a func- 
tion of Vg and B^. The relative phase of the two paths is 
^ — 2TT(j)/(/)Q, where cj) ~ B±Acs is the magnetic flux en- 
closed by the current flow and Acs = LcsWcs is the effective 
area. For our geometry, the latter depends mainly on the ef- 
fective width Wcs{Vg), which is a function of Vg through the 
energy dependence of A (fc^) (W^tf ^ W[coth{W/ X) - X/W] 
for X/W <C 1). As expected, the transmission is an simple 
oscillatory function of B±. Note that the shorter period cor- 
responds to the maximum effective area, 0o / (LcbW) ~ 1.5T 
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FIG. 3: (Color online) a) Density plot of the spin-resolved transmis- 
sion T++ as a function of the depth of the potential well Vg and the 
perpendicular magnetic field B± for a 32-ZGNR and ^isBa = t' /2, 
Ev = -QM, L = 2400a, A = 30a; b) same for T_ + ; c) magnetic 
field dependence of T++ (open symbols) and T_+ (filled symbols) 
for Vg/t = —0.53 (squares), —0.515 (circles) and —0.51 (triangles). 



of the edge states) there are no oscillations. The total trans- 
mission T — r_|-+ + T (. is constant, implying that the effect 

of the field is to produce a pure spin rotation. It is worth point- 
ing out that B± can not be too large to avoid a transition to a 
ferromagnetic state. 

Interestingly enough, there is also a way to produce a con- 
trolled spin rotation using an all-electrical-setup . The key is 
to change ^ by inducing a difference between the wavevectors 
of the two paths, and therefore changing their relative plane 
wave phase. This can be achieved by applying a small trans- 
verse electrical field that changes the energy of the two paths, 
and then the wavevectors, in a small fraction and in oppo- 
site directions — note that only a change Sk^ ~ 27r/i is re- 
quired. The transverse potential is described by adding a term 
Vt {[y^ - Wl2]2/W)f{x,) to Vgf{x,) in Eq. ©. Figure|4] 
shows the spin-dependent transmission for this setup. As for 
the previous case, there are clear oscillations indicating the 
rotation of the spin of the carriers, even for a very small trans- 
verse field Et = 2Vt/W (~ 2^lV/k for Vt ^ h IQ-H) . 
Again, the rotation disappears for Vg > V*. The period of 
the oscillations is in good agreement with the estimated value 

= (fc^ — k^)LcS where is the wave vector of the edge 
states with energy Ep — Vg + ri{VT) and (Vt) is the aver- 
age value of the transverse potential in the corresponding edge 
state {{Vt) ^ Vt[coI\v{N a / X) - X/Na] for X/W < 1). 

Adiabatic transport is not possible for ^isBa > t' as elec- 
trons reach a point where ~ before they penetrate the 
well and are then reflected. However, transport is still possi- 
ble due to a resonant mechanism that involves the upper 'w'- 
shaped edge states band (Fig. [U. This involves a Landau- 
Zener like transition between bands, so that the width of the 
resonances increases as the potential profiles is more abrupt. 
Some of those resonances are already apparent in Figs. [3]and 
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FIG. 4: (Color online) a) and b) Same as Fig. [3] but as a function 
of the transverse electrostatic potential Vr. ; c) transverse electric 

field dependence of r++ (open symbols) and T \. (filled symbols) 

for Vg/t = -0.544 (□,■), -0.53 (o,.) and -0.515 (A,A). 



m We note that in Fig. [3] they present a period of 2(/)o/Aofr. 
This is also present when the incident electron has its spin 
direction in z, where we would have naively expected no de- 
pendence with B± as the electrons in that case follow a single 
path. This is, however, not true as each minimum of the 'w'- 
shaped band involves edge states for electrons moving in one 
direction but extended states for those moving in the opposite 
direction — another unique characteristic of the ZGNR band 
structure. The phase difference is then related to half the area 
of the PW. A detailed analysis ll28ll shows that the spin rota- 
tion is still possible for some of the resonances, showing that 
the effect is robust against the precise value of Ba- 
in summary, we showed that ZGNRs present interesting in- 
terference phenomena in the presence of a PW. Moreover, the 
spin-dependent structure of the edge states allows for a con- 
trolled rotation of the spin of the carriers by either magnetic 
or electric fields. Since the characteristic of the zigzag termi- 
nation seems to be generic [31] and robust against disorder 
[l9L we expect these effects to manifest in less ideal samples, 
opening a new alternative for spintronics in graphene. 
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